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Abstract. We consider a generalized nonlinear Schrodinger equation (NLS) with a power nonlinearity |-i/>| 2 '*^/) of 
focusing type describing propagation on the ramified structure given by N edges connected at a vertex (a star graph). 
To model the interaction at the junction, it is there imposed a boundary condition analogous to the 8 potential of 
strength a on the line, including as a special case (a = 0) the free propagation. We show that nonlinear stationary 
states describing solitons sitting at the vertex exist both for attractive (a < 0, representing a potential well) and 
repulsive (a > 0, a potential barrier) interaction. In the case of sufficiently strong attractive interaction at the vertex 
and power nonlinearity n < 2, including the standard cubic case, we characterize the ground state as minimizer of 
a constrained action and we discuss its orbital stability. Finally we show that in the free case, for even N only, the 
stationary states can be used to construct traveling waves on the graph. 

PACS numbers: 05.45.Yv, 03.75.Kk, 42.65.Tg 



1. Introduction 

The Nonlinear Schrodinger Equation (NLS) is a well 
known model for several physical systems, and a para- 
digm for the behaviour of nonlinear dispersive equations 
PP. We recall two main fields of use of NLS in physics 
and its applications, to which we refer in the following: 
propagation of electromagnetic pulses in nonlinear media 
(typically laser beams in Kerr media or signal propaga- 
tion in optical fibers) , and dynamics of Bose-Einstein con- 
densates (BEC). It is well known that under some condi- 
tions (e.g. power nonlinearities) NLS admits symmetries, 
which in turn are associated to the existence of so called 
solitary solutions, usually in the form of standing or trav- 
elling waves. Symmetries are heavily affected and usu- 
ally broken by the presence of inhomogeneities of various 
type, but some of them persist and give rise to quite in- 
teresting phenomena, such as defect induced modes [3] , 
i.e. standing solutions strongly localized around the de- 
fect. For example the presence of defect modes affects 
propagation allowing trapping of wavepackets, as exper- 
imentally shown in the case of local photonic potentials 
in [3]; and on the other hand, nonlinearity can induce 
escaping of solitons from confining potentials, as demon- 
strated in [5], or enhanced quantum reflection of matter 
wave solitons [B]. A last interesting phenomenon is the 
strong alteration of tunneling through potential barriers 
in the presence of nonlinear defocusing optical media [7]. 
In this Letter we consider NLS propagation through junc- 
tions in networks showing analytically that defect or bet- 
ter "junction" modes arise, and giving their explicit con- 
struction and main properties. When the dynamics of a 
BEC takes place in essentially one dimensional substrates 
("cigar shaped" condensates) or a laser pulse propagates 
in optical fibers and thin waveguides, the question arises as 
to the effect of a ramified junction on propagation and on 
the possible generation of stable bound states. The study 



of the behavior of NLS on networks is a not yet devel- 
oped but growing subject. Concerning situations of direct 
physical interest we mention the analysis of scattering at 
Y junctions ("beam splitters") and other network configu- 
rations ("ring interferometers") for one dimensional Bose 
liquids discussed in [8]. Interplay between propagation 
in photonic T-junction and nonlinearity is studied in [§]. 
Some more results are known for the discrete chain NLS 
model (DNLS), see in connection with the present paper 
[TU] . Other recent developments are in [TT] and in |12) . 
where in particular stationary scattering from complex 
networks sustaining nonlinear Schrodinger propagation is 
studied. Again regarding scattering, we mention that the 
behaviour of fast solitons colliding with a Y-junction was 
studied in [13] . Here we would like to give a first rigorous 
determination of defect modes at junctions, i.e. nonlin- 
ear bound states, fully taking in account nonlinearity of 
the model NLS equation, and some properties of the cor- 
responding standing waves, in particular their stability. 
Even if the mathematical model of a ramified structure 
here used is oversimplified compared to a realistic exper- 
imental setup (but it is the common one e.g. in the the- 
oretical description of quantum wires), it allows to grasp 
the main effect we want to exhibit, i.e. the birth of gen- 
uinely nonlinear bound states. We recall that the linear 
Schrodinger equation on graphs has been for a long time 
a very developed subject [T3], due to its applications in 
quantum chemistry, nanotechnologies and more generally 
mesoscopic physics. We consider the simple configuration 
of a star graph Q made of N half lines with common origin. 
The associated Hilbert space for the Schrodinger dynamics 
is L 2 {Q) = (BjLi £ 2 (IR + )- Its elements will be represented 
as vectors functions with square integrable components, 
namely 
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To define a dynamics on the graph, we split it into a linear 
part described by a Hamiltonian operator and a nonlinear 
part. To assign the Hamiltonian operator on the graph 
one gives it on every edge; moreover a boundary condition 
at the vertex will be needed to assure a unitary dynamics. 
As an operator, the Hamiltonian will be correctly defined 
through both the action on every edge and the boundary 
condition at the vertex. So it will embody the parameters 
possibly present in the boundary condition. To be specific, 
we assume that the linear part of the Hamiltonian is the 
selfadjoint operator H a (see e.g. [Hj) which acts as the 
free Laplacian on every edge, (H a 1 i') i = —ip" , and with 
the following boundary conditions containing the param- 
eter a € R 



JV 



^i(O) 



\i>'h = o^i(o) 



(2) 



k = l 



H a generalizes to the graph the well known Schrodinger 
operator with delta potential of strength a on the line, 
and similarly to that case, the interaction is encoded in 
the boundary condition. For example the here mostly in- 
teresting case a < can be interpreted as the presence of 
a deep attractive potential well or attractive defect at the 
vertex, while a > represents a model of potential barrier 
at the junction. The case a — in ([2| plays a distinguished 
role and defines what is usually given the name of free (on 
the line one has continuity of wavefunction and its deriv- 
ative and interaction disappears) or Kirchhoff boundary 
condition. More general boundary conditions can be con- 
sidered to take into account properties of different physical 
models. To pose the NLS on the graph we define products 
componentwise, 



and Ifl 2 ^ 



(3) 



Now it is well defined the NLS equation on the graph, 

i|*(t) = H a *(t)-\*(t)\ 2 ^(t) (4) 

where /i > 0. This amounts to a system of N scalar NLS 
equations on the halfline, coupled through the boundary 
condition at the origin included in the domain ([2| . Notice 
that we consider a model with generic power nonlinear- 
ity; in most of physical applications the cubic (fx = 1) 
case is studied, but other nonlinearities appear in the lit- 
erature, in particular in the analysis of wave collapse in 
plasma physics and nonlinear optics when high nonlinear- 
ities (e.g. non Kerr) are present (see [H]). Concerning 
general properties of the equation, for /i > existence 
and uniqueness for initial data in the relevant state spaces 
holds ([131 [H]), an d if < \i < 2 solution exists for all 
times and no collapse occurs. Finally, as in the standard 
NLS on the line, mass M(*) = 1 1*111 and energy E[^f] are 



conserved, where 
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dx 



;iv>i(o)i 2 (5) 



2. Stationary states 



We are interested in the stationary solutions of the pre- 
vious equation, by which we mean solutions of Q of the 
form 

*(t,x) = e iut Hf u (x) . (6) 
The parameter uj is interpreted as propagation constant 
in nonlinear optics and as chemical potential in BEC. The 
amplitude *„ satisfies 

H a ^ u - |ig 2 ^ = -wtf w uj>0. (7) 

Locally H a acts as the Laplacian and so on every edge we 
must seek L 2 (R + ) solutions of the equation 

- 4>" - \4>\^ l (f) = -utf> u>0. (8) 

The most general one is given by the function 

i 

<j>(a, a; x) — a [(fi + 1)uj] 2fl sech" (fj,y/u(x — a)) (9) 

where |cr| = 1 and a 6 E. Therefore the components (* w )i 
of a stationary state * w must satisfy 

{y u ) i (x) = <f>(<7i,a i ;x). (10) 

Now we impose boundary conditions |2]). The continuity 
condition in ^ immediately implies a± = ... = fijv and 
a, = £id with Si = ±1 and a > 0. Due to the bell shape 
of the function <j), we shall say that in the i-th edge: there 
is a bump if ai > 0, that is, if Si = +1; there is a tail if 
a, < 0, that is, if Ej = —1. Now we determine e$ and a. 
The second boundary condition in ([2| implies 



N 

tanh(/z-\/wa) £j 

i=l 



a 

s/uj 



(11) 



Eq. (Ill gives as a first constraint that Yli=i £ i mus t have 
the same sign of a. That is for a > the stationary state 
must have strictly more bumps than tails while for a < 
we must have more tails than bumps. For every such a 



configuration condition (11) fixes uniquely a. We choose 



to index the stationary states by the number j of bumps. 
So *£ is given by 



0(a J ;x) i = l,...j 
(pi-a^ix) i = j + 1,...N 



: arctanh 



(2j - A0v^ 



(12) 



(13) 



In Eq. (12 I wc set <7j = 1 and omitted it in the argument 
of 4> to simplify notation. To summarize, solutions of ^ 
for a > are given by *£ with j = [N/2 + 1],...,N and 
for a < by with j = 0, . . . , [(N - l)/2] where [s] is 
the integer part of s. Stationary states for the three-edge 
star graph are portrayed in Fig. [TJ for the attractive case, 
and in Fig. [2j for the repulsive case. Notice that from 
( 13 1 follows the lower bound on existence of bound states 
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fpt < u) . We can interpret this fact for a < noticing 
that ([7]) neglecting nonlinearity is the eigenvalue equation 
for the linear part of the hamiltonian corresponding to en- 
ergy E = —lu, and taking in account the known fact that 

the linear graph hamiltonian H a has the ground state en- 
2 

ergy — j^; so the lower bound means that the nonlinear 
standing waves bifurcate from the vanishing wavefunction 
at the ground state energy of the linear problem. It can 
be rigorously proved that there are no solutions at all for 

3*>w • 





Figure 1. Nonlinear stationary states: 
a< 0,N = 3,j = 0, 1 





Figure 2. Nonlinear stationary states: 
a > 0,N = 3,j = 2,3 

Let us now consider the free case, a = 0. If N is odd the 
only solution to ( 11 ) is a = 0, and the stationary state 



is unique and it is given by 



This corresponds to N half solitons glued at the vertex. If 
N is even a family of solutions is given by 



a G 



tj>(-a,x), i = l,...N/2 
cj)(+a,x), i = N/2 + l,...N 

(15) 

Here a is a free parameter. In this case the graph can 
be considered just as a set of N/2 copies of the real line. 
The above solutions 4^ can be interpreted as N/2 iden- 
tical NLS solitary waves on each real line translated by a 
quantity a. We remark finally that from this last family 
of stationary states one can construct, through Galilean 
invariance and on each copy of the line, traveling waves of 
the form 



vt. 



(16) 



S U [*]=E[*] 



(17) 



i=l,...,N. 



(14) 



3. Ground state and variational properties 

Now we want to discuss existence and identification of the 
ground state of the system among the domain £, defined 
as the space of square intcgrablc functions with square in- 
tegrable derivative on every edge (sometimes indicated in 
the mathematical literature as H 1 (Q), the first Sobolev 
space on the graph) and continuous in the vertex. The 
stationary states of the NLS on graphs previously deter- 
mined are critical points (<S'^[ V I'] = 0) of the so called action 
functional S u 

LU 

2 

For a BEC the action is nothing but the grand poten- 
tial functional (at zero temperature) and the parameter 
lj, usually denoted as fi is the chemical potential. The 
ground state is the stationary state on which the action 
is minimal. For the main example of attracting defect 
this is the N tail state, but this does not mean that it is 
a minimum of the action. The action, due to the nega- 
tive nonlinear part of the energy, is unbounded from be- 
low, and to study its minima it is necessary to constrain 
it suitably. So we minimize S u on the so called natural 
constraint {* G £ s.t. 1^} = {S' u = 0} , which 
contains all the stationary states by definition, and where 
the action turns out to be bounded from below. One gets 
the following result, the non trivial proof of which will be 
given elsewhere (see [16] where the needed mathematical 
tools to get the result are given in detail). 

Theorem 1. Let [i > 0; there exists a* < such that for 

2 

a < a* and lu > , the action S u attains its constrained 
minimum on the N tail state \J r 2>. 

The limitation in a is due to the fact that for a = the 
iV-tails state is not a minimum but a saddle of the en- 
ergy functional. This is shown in |17j , constructing a trial 
wavefunction given by a quasi soliton almost completely 
concentrated on a single edge with an energy lower than 
the energy of a stationary state. To avoid this difficulty 
and guarantee the minimizing character of the stationary 
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state given in the Thcorcm[T]a sufficiently strong potential 
well is needed 

A second physically relevant constrained problem is to 
minimize the energy at constant mass, i.e. 

mf{E{q} s.t. * e £, M [*] = M[*°]} . (18) 

Here we only give an explicit example for the especially in- 
teresting case of the cubic NLS. A direct calculation shows 
that the mass of the bound states is independent of 
J : ll^illl = 2Nu> 2 + 2a . Notice that this degeneration is 
a special property of cubic case, not shared in general by 
different nonlinearities. This makes the energy spectrum 
at fixed mass exactly computable: 



EM 



N 



1 



-w 2 



3 (2j - TV) 2 



(19) 



The N bound states (formula holds true for both signs of 
a) have an energy increasing with the number of bumps. 
A closer analysis shows that this ordering holds true for 
the general power nonlinearity, but in this case it is not 
possible to give an explicit nonlinear energy spectrum. 
We end this section noticing that for a < the only sta- 
tionary state respecting the symmetry of the interaction at 
the vertex is the ground state, while in the excited states 
the symmetry with respect to the edges is broken. 



4. Stability of the ground state 

Now we come to the stability of the ground state. Due to 
the U(l), or phase, invariance of Q the stability has to 
be considered as orbital stability [21 [T] [TB] , which means 
ordinary Lyapunov stability up to symmetries: a solution 
remains as close as desired to the orbit {e lujt of a 
ground state if it starts close enough to it. Establishing 
orbital stability is generally more difficult of establishing 
stability of an isolated equilibrium point. An analogous 
study for the delta potential on the line was given in [2], 
second reference. As in the scalar case, the NLS on a 
graph turns out to be a hamiltonian system on the real 
Hilbert space of the couples of real and imaginary part of 
the wavefunction. Decomposing = u + iv, Q is equiv- 
alent to the canonical system 



d 
dt 



JE'[u,v] 



J 



I 

-I 



(20) 



where the Hamiltonian E[u, v] is obtained from ^ after 
posing = u + iv . Linearization of the hamiltonian sys- 



tem (201 around the stationary state is achieved by sub- 
stituting (^/(t))k — k + i]k + ipk)e wt and neglecting 
higher order terms than linear in ( 20 1 . The real vector 
functions t] and p satisfy 



d 
dt 



(21) 



where L — diag(£_,£+) and 
d 2 



(£-),.* = - 



_ + w -|<*l a/i K* 



d?_ 

dx- 



+ u,-(2 M + l)|< fc | 2 ''W 



(22) 



£_ and £ + are matrix self adjoint operators acting on the 
real vector functions 77 and p satisfying ([2]). Note that 
= 0. Precise conditions to have orbital stability 
for general hamiltonian systems and in particular for NLS 
equations are now classical [HI [19] . They are based on the 
properties of linearization C and on the so called Vakhitov- 
Kolokolov (|20j) criterion. These results imply in partic- 



ular that solitary solutions of (20) are orbitally stable if 
i) spectral conditions hold: i\) ker£ + = {W"} and the 
rest of the spectrum is positive; £2) n(C-) — 1 where the 
left hand side is the number of negative eigenvalues, ii) 



Vakhitov-Kolokolov condition ^?||^| 



I > holds. 



Theorem 2. Let p, € [0, 2], a < a* < 0, cj > 
the ground state \t„ is orbitally stable in £. 



N 2 ■ 



The 



Concerning Vakhitov-Kolokolov condition one has 
d 



dui 



+ 



= C 



f 1 (I-; 2 );- 1 dt 







[G- 





(23) 



1 - 



N 2 o 



with C = C{N, p, u) 



> 



Concerning ii), 



'JjJJ is positive, and we prove now that it is the (simple) 
ground state of C+. An integration by parts allows to 
rewrite 



(£+*,*) 



E 

k=l 



+00 . 



01 



dx y *°. 



! da; (24) 



where the finite term in the integration is vanishing due 
to the S boundary condition. So (£+<!', "J) > for every 
\& not coinciding with ^j^, which is the only eigenvector 
with eigenvalue 0. Concerning 12), note that the ground 
state is a strict local minimum of S u constrained on the 
codimension one natural manifold 7 L J 1 (0). So C is positive 
on the tangent space to Z~ 1 (0) at the ground state. Corre- 
spondingly, C could have at most one negative eigenvalue. 
In fact it has one, being , *°) < 0. This ends the 

proof. 

We add some remarks. In the first place notice that from 



formula ( 23 ) it follows that for p > 2 there exists u* such 
that is orbitally stable for ui G an d is or- 

bitally unstable for u> > ui* .So one has orbital stability 
of the branch of ground states also in the supercritical 
regime p > 2, where global existence of solutions is not 
guaranteed and the solution could blow up in finite time, 
if the frequency uj of the standing wave is not too high. 
Nothing rigorous is known about orbital stability or insta- 
bility of excited states; it is expected that excited states 
are unstable. 
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5. Conclusions 

We begin with this letter the analysis of nonlinear wave 
propagation on ramified structures and networks. We have 
considered a star graph (N halflines, "edges", attached 
at a single vertex) which sustains a nonlinear focusing 
Schrodinger dynamics of power type. Interaction between 
the fields on different edges is through a boundary con- 
dition at the vertex. For this system explicit nonlinear 
bound states exist when certain classes of boundary con- 
ditions (5 type, attractive or repulsive) are considered. In 
the most common case of Kirchhoff vertex, for an even 
number of edges only, traveling waves exist. Let us stress 
the fact that there exist excited states beside the ground 
state, For example they give rise to a nonlinear energy 
spectrum for the system, here displayed in the case of cu- 
bic nonlinearity. Moreover, the ground states have a vari- 
ational characterization and they are orbitally stable for 
sufficiently attractive vertex and nonlinearity below a crit- 
ical value. The study can be extended to the NLS with dif- 
ferent nonlinearities, for example to combination of several 
power nonlinearities, such as the important case of cubic- 
quintic case. Moreover similar techniques can be applied 
to other nonlinear dispersive equations admitting solitons 
solutions or standing waves. Finally, apart from the pos- 
sible applications of the model to the description of real 
systems involving propagation on networks, the obtained 
results rise several general interesting issues. We mention 
two of them, which seem to be important for a general 
comprehension of the dynamics of solitary waves. The 
first is the analysis of stability and instability of excited 
states of the system. It is known that excited standing 
waves exist in many nonlinear models, but it is in general 
quite difficult to study them due to the rather implicit 
information at disposal. The present model has explicit 
excited states, and it appears as promising their stability 
study, in the first place at a linearized level (given by Eq. 
(21)) , and then concerning orbital (in)stability. A sec- 
ond interesting issue concerns the asymptotic stability of 
solitons, i.e. the fact (to be proven in every single model) 
that the evolution of an initial state in a suitable class 
asymptotically decomposes in one or more solitary waves, 



radiation driven by the linear evolution and a small re- 
mainder. When true this property is obviously quite im- 
portant from the physical point of view because it yields 
to a strong simplification of the behaviour of the system 
in the asymptotic regime. A possible proof is related in 
the first place to the spectral properties of linearization 
J d and to the decay properties of solutions of the lin- 



earized equation (21), and again one can hope to exploit 



the rather explicit expressions of these objects. 
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